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ABSTRACT 

"in  this  paper  we  survey  recent  results  concerning  global  existence  and 
decay  of  smooth  solutions  of  certain  quasilinear  hyperbolic  Volterra  equations 
which  provide  models  for  the  motion  of  one-dimensional  viscoelastic  solid  of 
the  Boltmsann  tyr>e .  We  also  sketch  the  derivation  of  these  equations  from 
physical  principles,  discuss  the  physically  appropriate  assumptions ,  and  prove 
a  special  case  of  a  new  existence  theorem  for  the  Cauchy  problem. * 
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GLOBAL  EXISTENCE  AND  ASYMPTOTICS  IN 
ONE-DIMENSIONAL  NONLINEAR  VISCOELASTICITY 


W.  J.  Hrusa 1 ' 2  and  J.  A.  Nohel1 

1.  Introduction 

For  nonlinear  elastic  bodies,  the  balance  laws  of  continuum  mechanics 
lead  to  equations  of  notion  (of  hyperbolic  type)  which  have  the  property  that 
smooth  solutions  may  break  down  in  finite  time  due  to  the  formation  of  shock 
waves •  Some  material  models  of  physical  interest  incorporate  a  nonlinear 
"elastic-type"  response  in  conjunction  with  a  natural  dissipative  mechanism. 
For  such  materials  it  is  important  to  understand  the  effects  of  dissipation  on 
solutions  of  the  equations  of  motion. 

Some  dissipative  mechanisms  (e.g. ,  viscosity  of  the  rate  type  in  one 
space  dimension)  are  so  powerful  that  globally  defined  smooth  solutions  exist, 
even  for  very  large  initial  data.  A  much  more  subtle  type  of  dissipation,  due 
to  memory  effects,  arises  in  viscoelasticity  of  the  Boltmann  type. 

In  this  paper  we  discuss  global  existence  and  decay  of  smooth  solutions 
of  certain  quasllinear  hyperbolic  Volterra  equations  which  provide  models  for 
the  motion  of  nonlinear  viscoelastic  solids  of  the  Boltzmann  type.  In  Section 
2  we  formulate  the  dynamic  problems  to  be  considered  and  discuss  the  relevant 
assumptions.  In  Section  3  we  give  a  survey  of  known  results.  Theorem  3.1  is 
newr  the  complete  proof  will  appear  elsewhere  [17].  Finally,  in  Section  4,  we 
prove  a  special  case  of  Theorem  3.1. 
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«•  restrict  our  Attention  throughout  to  one-diaensional  notions* 

Although  the  details  have  not  been  carried  out  completely,  analogous  results 
can  be  obtained  for  nultidinensional  viscoelastic  solids  of  the  Boltzmann 
type*  This  is  discussed  briefly  in  (?].  Local  existence  results  which  are 
applicable  to  multidimensional  bodies  occupying  all  of  space  have  been  given 
by  Grimmer  and  Zeman  (12]. 

He  close  the  introduction  with  some  remarks  on  notation.  Let  D  be  a 
subset  of  R  *  R.  For  a  function  w  «  D  ♦  R  we  use  subscripts  x  and  t 
(or  t)  to  indicate  partial  differentiation  with  respect  to  the  first  and 
second  argument,  respectively.  Moreover,  we  use  the  same  symbol  w  to  denote 
the  mapping  t  **  w(*,t)  when  there  is  no  danger  of  confusion.  A  prime  is 
used  to  denote  the  derivative  of  a  function  of  a  single  variable,  and  the 
symbol  :•  indicates  an  equality  in  which  the  left  hand  side  is  defined  by 
the  right  hand  side.  All  derivatives  should  be  interpreted  in  the  sense  of 
distributions. 


2.  Formulation  of  Dynamic  Problems 

Consider  the  longitudinal  motion  of  a  homogeneous  one-dimensional  body 
that  occupies  the  interval  B  in  a  reference  configuration  (which  we  assume 
to  be  a  natural  state)  and  has  unit  reference  density.  We  denote  by  u(x,t) 
the  displacement  at  time  t  of  the  particle  with  reference  position  x 
(i.e.,  x  ♦  u(x,t)  is  the  position  at  time  t  of  the  particle  with  reference 
position  x),  in  which  case  the  strain  is  given  by 

e(x,t)  :■  ux(x,t)  .  (2.1) 

For  smooth  displacements,  the  equation  of  balance  of  linear  momentum  here 
takes  the  form 

utfc(x,t)  *  Ox(x,t)  +  f(x,t),  x  e  B,  t  >  0  ,  (2.2) 
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where  o  is  the  stress  and  f  is  the  (known)  body  force.  Equation  (2.2) 
must  be  supplemented  with  a  constitutive  assumption  (stress-strain  relation) 
which  characterizes  the  type  of  material  composing  the  body. 

If  the  body  is  elastic,  then  the  stress  depends  on  the  strain  through  a 
constitutive  equation  of  the  form 

0(x,t)  -  +(e(x,t))  ,  (2.3) 

where  4  is  an  assigned  smooth  function  with  +(0)  -  0,  and  the  resulting 
equation  of  motion  is 

utt  “  ♦‘Vx  +  f  ‘  <2*4> 

Experience  indicates  that  stress  increases  with  strain,  at  least  near 
equilibrium,  so  it  is  natural  to  assume  that  +' (0)  >  0.  Lax  [18]  and  MacCamy 
and  Mizel  [23]  have  shown  that  (2.4)  (with  f  =  0 )  does  not  generally  have 
globally  defined  smooth  solutions  no  matter  how  smooth  (and  small)  the  initial 
data  are. 

For  viscoelastic  materials  of  the  rate  type,  the  stress  depends  on  the 
strain  rate  as  well  as  the  strain.  A  simple  model  corresponds  to  the 
constitutive  relation 

o(x,t)  *  +(e(x,t))  +  Xet(x,t)  ,  (2.5) 

where  $  is  as  above  and  X  is  a  positive  constant,  which  leads  to  the 
equation 

»tt  ■  ♦'Vx  *  Xu,tx  *  '  •  <2-‘» 

Greenberg,  MacCamy,  and  Mizel  [11]  have  shown  that  the  Dirichlet  initial¬ 
boundary  value  problem  for  (2.6)  has  a  unique  globally  defined  smooth  solution 
provided  that  the  initial  data  are  sufficiently  smooth.  Viscosity  of  the  rate 
type  is  so  powerful  that  global  smooth  solutions  exist  even  if  the  initial 
data  are  very  large.  Similar  results  for  more  general  viscoelastic  materials 
of  the  rate  type  have  been  obtained  by  Dafermos  [3]  and  MacCamy  [19]. 
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Experience  indicates  that  in  certain  materials,  the  etreee  at  a  material 
point  x  depends  on  the  entire  temporal  history  of  the  strain  at  x.  Zn 
1876,  Boltzmann  [1]  proposed  the  constitute  ..on 

C(x,t)  -  cc(x, t)  -  m(s)C(x,t-s)ds  ,  (2.7) 

where  c  is  a  positive  constant  and  m  is  positive,  decreasing,  integrable, 
and  satisfies 

c  -  m(s)ds  >  0  .  (2.8) 

The  history  of  the  strain  up  to  time  t  ■  0  is  assumed  to  be  known. 

The  constant  c  measures  the  instantaneous  response  of  stress  to  strain, 
and  the  first  two  conditions  on  m  say  that  the  stress  "relaxes"  as  time 
increases  and  that  deformations  Which  occurred  in  the  distant  past  have  less 
influence  on  the  present  stress  than  those  which  occurred  in  the  recent 
past.  Equation  (2.8)  also  has  an  important  mechanistic  interpretation.  In 
statics,  i.e.  o(x,t)  =  o(x)  and  e(x,t)  =  e(x),  equation  (2.7)  reduces  to 

o(x)  ■  (c  -  Jq  m(s)d s)e(x)  ,  (2.9) 

and  thus  (2.8)  states  that  the  equilibrium  stress  modulus  is  positive. 

A  natural  nonlinear  generalization  of  (2.7)  is  provided  by  the 


constitutive  equation 

o(x,t)  «  +(e(x,t))  -  Jq  m(s)t(c(x,t-s))ds  (2.10) 

where  6  end  6  are  assigned  smooth  functions  with 

♦(0)  -  p(0)  -  0,  ♦'(0)  >  0,  t'(0)  >  0  ,  (2.11) 

and  m  is  positive,  decreasing,  integrable,  and  satisfies 

♦'(0)  -  (Jq  m(s)ds)p' (0)  >  0  .  (2.12) 

It  is  convenient  to  define  the  relaxation  function  a  by 

»(t)  »*  J*  m(s)ds,  t  8  (0,-)  ,  (2.13) 

and  the  equilibrium  stress  function  x  by 

X(C)  *«  6(C)  -  a(0)p(C),  (II  .  (2.14) 
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If  a  satisfies  the  preceding  conditions  then  a  is  positive,  decreasing, 
and  convex,  and  x'(0)  >  0. 

We  note  that  a*  =  -a.  Thus  (2.10)  can  be  written  in  the  fora 

o(x,t)  ■  4(c(x,t))  ♦  Jq  a'(s)+(e(x,t-s))ds  (2.15) 

and  (letting  t  »»  t-s)  also  in  the  fora 

o(x,t)  -  +(c(x,t))  +  J ^  a* (t-T)t(e(x,t))dr  .  (2.16) 

Ihe  corresponding  equation  of  aotion  is 

utt(x#t)  “  *tux(x,t»x  +  a,(t-r)t(ux(x,T))jcdT  (2.17) 

+  f (x, t ) ,  x  e  8,  t  >  0  . 

Observe  that  a',  rather  than  a,  appears  in  equations  (2.15),  (2.16), 
and  (2.17).  In  this  paper,  we  are  nomalizing  a  so  that  a(t)  ♦  0  as 
t  ♦  «•.  (See  (2.13).)  The  reader  is  cautioned  that  other  noraali nations  are 


frequently  used. 

An  appropriate  dynamic  problem  is  to  determine  a  smooth  function 
u  i  8  *  (-•»•)  *  *  which  satisfies  equation  (2.17)  lor  t  >  0,  together 
with  suitable  boundary  conditions  if  B  is  bounded,  and 

u(x,t)  •  v(x,t)  x  e  B,  t  <  0  ,  (2.18) 

where  v  is  a  given  smooth  function.  The  history  value  problem  (2.17), 

(2.18)  can  be  reduced  to  an  initial  value  problem  as  follows.  Define  a  new 
forcing  function  g  by 

g(x,t)  i»  f ( x, t )  +  a' (t-T)t(vx(x,t))xdT,  (2.19) 

x  e  B,  t  <  0  , 


and  initial  data  u0,  u1  by 

u0(x)  ■  v(x,0),  Uj(x)  “  vt(x,0),  x  e  B  .  (2.20) 

It  is  clear  that  u  is  a  solution  of  (2.17),  (2.18)  if  and  only  if  it  is  a 


solution  of  the  initial  value  problem 


(2.21) 


Utt(x,t)  -  ♦(Ux<X,t))x  +  J*  *’(t-T)f(ttx(X,T))xdT 

+  g(x,t),  x  6  B,  t  >  0  , 

u(x,0 )  -  Uq(x),  ut(x,0)  »  u^x),  x  e  B  .  (2.22) 

Conversely,  the  initial  valua  problan  (2.21),  (2.22)  can  ba  converted  to 
a  history  valua  problem  of  tha  form  (2.17),  (2.18)  by  constructing  suitable 
functions  v  and  f.  (Of  course,  such  a  procedure  does  not  uniquely 
determine  a  history  value  problem.)  For  consistency,  we  state  all  results  for 
initial  value  problems.  Clearly,  there  are  analogous  statements  for  history 
value  problems. 

We  consider  pure  initial  value  problems  (Cauchy  problems)  with  B  ■  K, 
as  well  as  initial-boundary  value  problems  with  B  ■  [0,1]  and  boundary 
conditions  of  Dirichlet,  Neumann,  or  mixed  type,  i.e. 

u(0,t)  -  u( 1 , t )  -  0,  t  >  0  ,  (2.23) 

ux(0,t)  •  ux(1,t)  -  0,  t  >  0  ,  (2.24) 

or 

u(0,t)  -  ux(1,t)  -  0,  t  >  0  .  (2.25) 

The  physical  interpretation  of  (2.23)  is  clear.  Under  certain  appropriate 
conditions,  (2.24)  is  equivalent  to 

0(0, t)  -  0(1, t)  -  0,  t  >  0  .  (2.26) 

See,  for  example,  [7].  (A  similar  comment  applies  to  (2.25).) 

For  initial-boundary  value  problems,  the  initial  data  and  g  should  be 
compatible  with  the  boundary  conditions.  For  example,  suppose  that  u  is  a 
classical  solution  of  (2.21),  (2.22),  (2.23),  (with  g  S  0  for  simplicity)  on 
[0,1]  x  [0,T]  for  some  T  >  0.  Differentiating  (2.23)  twice  with  respect 
to  t  yields 

ut(0,t)  -  ut(1,t)  -  utt(0,t)  -  utt(1,t>  -  0  (2.27) 


V  t  e  [0,T] 


If  (2.21),  (2.22),  (2.23),  and  (2.27)  ara  to  hold  at  t  -  0  (and  +'(u^) 


does  not  vanish),  then  Uq  and  u^  must  satisfy 

u0(0)  -  u0(1)  -  u1 (0 )  -  Ul(1)  -  uj(0)  -  u«(1)  -  0  .  (2.28) 

Violation  of  the  above  condition  should  be  interpreted  as  a  singularity  in  the 
initial  data  on  the  boundary.  Due  to  the  hyperbolic  nature  of  equation 
(2.21),  such  a  singularity  would  try  to  propagate  away  from  the  boundary  and 
into  the  interior.  Analogous  compatibility  conditions  are  required  for  (2.24) 
and  (2.25).  (If  g  f  0,  then  the  compatibility  conditions  also  involve  g.) 

3.  Survey  of  Results 

observe  that  if  a*  vanishes  identically,  then  (2.21)  reduces  to  an 
undamped  quae il inear  wave  equation.  If  a*  t  0  and  the  appropriate  sign 
conditions  are  satisfied,  the  memory  term  in  (2.21)  induces  a  weak  type  of 
dissipation.  A  great  deal  of  information  concerning  the  strength  of  this 
dissipative  mechanism  is  contained  in  the  work  of  Coleman  and  Gurtin  [2]  on 
the  growth  and  decay  of  acceleration  waves  in  materials  with  memory.  Roughly 
speaking,  they  showed  that  (under  physically  natural  assumptions)  the  ampli¬ 
tude  of  a  certain  type  of  weak  singularity  (involving  jump  discontinuities  in 
second  derivatives  of  u)  decays  to  zero  as  t  ♦  •,  provided  its  initial 
amplitude  is  sufficiently  small.  On  the  other  hand,  the  amplitude  of  such  a 
singularity  may  become  infinite  in  finite  time  if  its  initial  amplitude  is  too 
large. 

This  suggests  that  (2.21)  should  have  globally  defined  smooth  solutions 
for  sufficiently  smooth  and  small  data,  and  that  smooth  solutions  can  develop 
singularities  in  finite  time  if  the  data  are  suitably  large.  (Here  we  use  the 
term  data  to  mean  initial  data  and  forcing  function.)  Results  of  this  type 
have  been  obtained  by  a  number  of  authors.  At  the  present  time,  the  situation 
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concerning  existence  of  global  solutions  for  small  data  is  quite  well 
understood;  less  is  known  about  the  formation  of  singularities.  It  should  be 
noted  that  several  important  ideas  used  in  the  analysis  of  (2.21)  were 
motivated  by  the  work  of  Nishida  [27]  and  Matsumura  [26]  on  quasilinear  wave 
equations  with  frictional  damping. 

Local  existence  of  smooth  solutions  to  (2.21)  can  be  established  by  more 
or  less  routine  procedures.  (See.  for  example.  [7].)  The  local  arguments 
require  only  positivity  of  +'  and  smoothness  of  ♦,  ♦»  a,  and  the  data.  In 
particular,  they  are  insensitive  to  the  "sign"  of  the  memory  term  and  the  size 
of  the  data.  However,  rather  delicate  a  priori  estimates  are  needed  to  show 
that  local  solutions  can  be  continued  globally.  These  estimates  rely 
crucially  on  the  memory  term  having  the  correct  sign  and  the  data  being  small. 

For  the  special  case  f  =  global  existence  theorems  have  been 
established  by  MacCamy  [21].  Dafermos  and  Nohel  [6],  and  Staffans  [30].  In 


order  to  simplify  our  discussion  of  these  results,  let  us 
assume  that  g  2  0  and  consider  the  problem 

utt<X,t)  “  *(ux(x,t))x  +  JS  *,(t"T)^(ux(x'T))xdT  •  (3.1) 

x  e  B,  t  >  0  , 

u(x,0)  »  Uq ( x) ,  ut(x,0)  -  u ^ ( x ) ,  x  e  B  .  (3.2) 

The  main  hypotheses  on  ♦  and  a  are 

♦  e  C3(»),  <M0)  -  0,  ♦MO)  >  0  ,  (3.3) 

a,  a',  a"  e  L^O,-)  ,  (3.4) 

'  a  is  strongly  positive  definite  ,  (3.5) 

a(0 )  <  1  .  (3.6) 


(Some  additional  technical  assumptions  on  a  are  used  in  [21]  and  [6].)  We 
refer  the  reader  to  [29]  and  [30]  for  properties  of  strongly  positive  definite 
kernels.  We  note,  however,  that  twice  continuously  differentiable  a  which 


satisfy 


(-1>ka<)c)(t>  >0  V  t  >  0,  k  -  0,1, 2}  a*  $  0  ,  (3.7) 

are  automatically  strongly  positive  definite.  (Corollary  2.2  of  [29].) 
Condition  (3.6),  together  with  +*(0)  >  0,  simply  states  that  x'(°)  >  0. 

Remark  3. 1 :  We  note  that  a*  rather  than  a  appears  in  equation  (3.1).  Our 
normalizations  of  a  (with  a(*>)  -  0)  is  different  from  that  used  in  [21], 
[6],  and  [30].  For  this  reason,  the  conditions  on  a  above  are  in  a  slightly 
different  form  than  in  [21],  [6],  and  [30]. 


The  assumptions  needed  on  u0  and  u^  vary  slightly  depending  on  the 
type  of  boundary  conditions.  Roughly  speaking  it  is  required  that 

Uq,  u£»  “g***  uJ,Uj  €  L2(B>  (3.8) 

2 

and  that  the  L  (B)  norms  of  the  functions  listed  in  (3.8)  be  sufficiently 


small.  In  addition,  the  data  must  be  compatible  with  the  boundary  conditions 

2 

if  B  is  bounded.  It  is  not  assumed  that  the  L(B)  norm  of  uQ  is 
small.  However,  for  certain  initial-boundary  value  problems,  this  is  implied 
by  the  Po In car 6  inequality  and  smallness  of  the  L(B)  norm  of  u^. 


Under  the  above  assumptions,  the  initial  value  problem  (3.1),  (3.2), 

2 

with  B  *  R,  has  a  unique  solution  u  8  C  (R  *  [0,**))  such  that 


ut 'UX ,utt 'Utx 'uxx 'uttt 'uttx 'Utxx 'uxxx  e  CU0,-)|  L  <*> )  .  (3.9) 


Moreover,  as  t  ♦  •, 

Utt'  Utx'  uxx  +  0  in  1,2 lU)  '  0.10) 

u,,  u  ,  u^.  ,  u^  ,  u  ♦  0  uniformly  on  R  .  (3.11) 

t  x  tt  tx  xx 

Similar  conclusions  hold  for  initial-boundary  value  problems  for  (3.1)  with 
B  ■  [0,1]  and  boundary  conditions  (2.23),  (2.24),  or  (2.25).  The  precise 


decay  statement  depends  on  the  boundary  conditions.  For  (2.23)  or  (2.25) 
(i.e.,  Dirichlet  or  mixed  conditions). 


(3.12) 


u,  ut,  ux.  utt,  utx,  Uxx  ♦  0  uniformly  on  (0,1] 
as  t  ♦  while  for  (2.24)  (Neumann  conditions), 

V  utt,  utx,  uxx  ♦  0  uniformly  on  (0,1]  (3.13) 

as  t  ♦  •.  The  difference  is  due  to  the  fact  that  nontrivial  rigid  motions 
are  possible  under  (2.24),  but  not  under  (2.23)  or  (2.25).  See  (21],  (6],  and 
[30]  for  the  proofs.  (The  boundary  conditions  (2.25)  are  not  discussed 
explicitly,  but  the  same  proofs  apply  with  only  trivial  modifications.) 

Remark  3.2:  If,  under  boundary  conditions  (2.24),  it  is  assumed  that  the  data 
have  zero  average  spatially  then  the  solution  will  have  zero  average  spatially 
and  (3.13)  can  be  replaced  by  (3.12).  k  Neumann  problem  can  always  be  reduced 
to  one  in  rtiich  the  data  have  zero  average  by  superposition  of  a  rigid 
motion.  (See,  for  example,  [7]  or  [16].) 

Remark  3.3t  The  above  results  remain  valid  if  a  suitably  smooth  and  small 
forcing  function  g  (which  behaves  properly  as  t  ♦  •)  is  included  in 

(3.1) .  See  [6],  [21],  and  [30].  (See  also  Theorem  3.1  below  for  an 
indication  of  the  type  of  assumptions  required  of  g. ) 

On  the  other  hand,  Hattori  [13]  has  shown  that  if  ♦’(C)  >  0 
V  £  e  R  and  4"  ?  0#  then  there  are  smooth  initial  data  (compatible  with  the 
boundary  conditions)  for  which  the  initial-boundary  value  problem  (3.1), 

(3.2) ,  (2.23),  with  B  -  [0,1],  does  not  have  a  globally  defined  smooth 
solution.  Such  data  must  necessarily  be  large  in  view  of  the  aforementioned 
existence  results.  The  precise  manner  in  which  loss  of  regularity  occurs  is 
not  discussed  in  [13].  Markowich  and  Renardy  [25]  have  obtained  numerical 
evidence  which  indicates  the  formation  of  shock  fronts  in  smooth  solutions  of 
the  initial  value  problem  (3.1),  (3.2)  with  B  *  R  and  suitably  large  Initial 
data. 
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The  following  idea  of  MacCamy  reveals  that  there  is  a  close  similarity 


between  (3.1)  and  a  wave  equation  with  frictional  damping.  Observe  that 
$(  u  )  can  be  expressed  in  terms  of  u**  through  an  inverse  linear  Volterra 
operator.  An  integration  by  parts  can  then  be  used  to  transfer  a  time 
derivative  from  to  the  resolvent  kernel  associated  with  a*.  This 

introduces  a  frictional  damping  term  and  renders  the  memory  term  a  linear 
perturbation  of  lower  order. 

More  precisely,  the  (scalar)  linear  Volterra  operator  L  defined  by 

( Lw ) ( t )  w(t)  +  jjj  a' (t-T)w(T)dr,  t  >  0  ,  (3.14) 

is  invertible  with  inverse  given  by 

(L_1w)(t)  -  w(t)  +  J*  k(t-T)w(T)dt,  t  >  0  ,  (3.15) 

where  k  is  the  resolvent  kernel  associated  with  a',  i.e.  k  is  the  unique 

solution  of 

k(t)  +  Jq  a' (t-T )k(T)dx  *  -a'(t),  t  >  0  .  (3.16) 

Using  (3.15)  to  solve  (3.1)  for  ♦<ux>x  in  terms  of  ufct  yields 

$<ux(x,t))x”  utfc(x,t)  +  Jq  k(t-T)utt(x,x)dr  (3.17) 

x  e  B,  t  >  0  . 

After  an  integration  by  parts,  this  becomes 

utt(x,t)  +  k(0)ut(x,0)  *  *<ux(x,t))x  +  k(t)u1(x)  (3.18) 

-  Jp  k' (t“T )ut(x,T)dT  ,  X  e  B,  t  >  0  , 

where  use  has  been  made  of  (3.2).  It  follows  from  (3.16)  that  k(0)  « 

-a*(0),  and  thus  the  term  k(0)ut  has  a  damping  effect  if  a*(0)  <  0. 

This  form  of  the  equation  is  extremely  convenient  for  many  purposes. 

Remark  3.4:  If  f  =  ♦,  then  (2.21)  also  arises  in  a  mathematical  model  for 
heat  flow  in  materials  with  memory.  For  the  heat  flow  problem,  (3.3),  (3.4), 
and  (3.5)  are  still  appropriate,  but  (3.6)  should  be  replaced  by  a(0)  «■  1. 
This  seemingly  minor  change  leads  to  major  differences  in  the  analysis.  The 


-11- 


.  C|  ■  ISMWJ.V.WWUWS  -  -  ---  -—>■  "“-'I  ■  ...  ...  ..  .**1 


memory  term  actually  has  a  slightly  stronger  dissipative  effect  in  this 
situation.  (See  (20],  [6],  and  [30].) 


For  the  general  case  with  4  different  from  4,  Dafermos  and  Mohel  [7] 
exploited  the  positivity  of  x'(°)  and  the  strong  positive  definiteness  of 
a  to  obtain  global  a  priori  estimates  for  solutions  of  initial-boundary  value 
problems  with  B  *  [0,1].  They  integrate  by  parts  and  use  (2.14)  and  (2.22) 
to  rewrite  (2.21)  in  the  form 

u  (x,t)  -  x<u  (x#t))„  ♦  J«  a(t-T )#(u  (x,T) )  dr  (3.19) 

+  a(t)4' (uq(x) )uj(x)  +  g(x,t)  , 
x  e  b,  t  >  o  . 

They  obtain  estimates  for  certain  higher  order  derivatives  directly  from 
(3.19)  and  use  the  Poincari  inequality  to  estimate  lower  order  derivatives. 
Their  procedure  yields  global  existence  (and  decay)  of  smooth  solutions  for 
small  data  with  B  »  [0,1]  under  boundary  conditions  (2.23),  (2.24),  or 
(2.25).  However,  due  to  the  lack  of  PoincarS-type  inequalities  on  all  of 
space,  their  results  do  not  apply  to  the  pure  initial  value  problem  (2.21), 
(2.22)  with  B  -  *. 

Regarding  4,  4#  and  a,  they  assume  that 

4.  4  «  C3(«),  4(0)  -  4(0)  -  0  ,  (3.20) 

4' (0 )  >  0,  4* (0)  >  0,  x'(°)  >  0  »  (3.21) 

and  that  (3.4)  and  (3.5)  hold.  Their  assumptions  on  the  data  and  the 
conclusions  of  their  existence  theorems  are  essentially  the  same  as  those 
stated  previously  for  initial-boundary  value  problems  in  the  special  case 
4  3  4. 

Subsequently,  Hrusa  and  Nohel  [17]  established  a  global  existence 
theorem  for  the  Cauchy  problem 
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(3.22) 


utt<x,t)  -  ♦(ux(x,t))x  +  J*  a,(t-T)t(ux(x#T))x<lT 
+  g(x,t) ,  x  e  ft,  t  >  0  , 

u(x,0)  ■  uQ (x) ,  ut(x,0)  -  u1(x),  x  e  R  .  (3.23) 

He  state  a  slightly  simplified  version  of  this  result. 

Theorem  3.1;  Assume  that  (3.20)/  (3.21),  (3.4),  (3.5)  hold,  and  that 
a  satisfies  some  (mild)  additional  technical  conditions,  Then,  there  exists 
a  constant  u  >  0  such  that  for  each  Ug,  u^  :  R  ♦  *  and  g  :  R  x  [0,«)  ♦  R 
with 


“o  e  "O'  "o'  "}”'  V  «i'  "1  «  I'2<“  ' 

(3.24) 

g,  gfc,  gx  e  c([0,«)»  l2 (r)  )  , 

(3.25) 

g,  gt  e  L1  ( (0,*)  |  r,2(*))  , 

(3.26) 

9x'  gtt  e  ' 

(3.27) 

oil.  tug ( ic)2  +  Ug(x)2  +  u-'(x)2  ♦  u1(x)2 

(3.28) 

+  u{(x)2  +  u"(x)2}(x)dx)1^ 

+  sup  (J*  {g2  +  g2  +  g2}(x,t)dx)^ 

t>0  r  X 

+  j"  (J^  (g2  ♦  g2}(x,t)dx)^  dt 

+  < Jq  C  +  9^t^x*t)dxdt)/^ 

<  M  , 

the  initial  value  problem  (3.22),  (3.23)  has  a  unique  solution 
u  6  C2(R  x  (0,«) )  which  satisfies  (3.9).  Moreover,  as  t  ♦  •,  (3.10)  and 

(3.11)  hold. 

The  proof  combines  certain  estimates  of  Dafermos  and  Hohel  (7]  for  higher 
order  derivatives  (which  remain  valid  for  B  -  ft)  with  a  variant  of  MacCamy's 
procedure.  (See  [17]  for  the  details.)  The  additional  technical 
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assumptions  on  a  (which  are  stated  precisely  in  (17])  are  not  very 
restrictive;  their  purpose  is  to  ensure  integrability  of  certain  resolvent 
kernels.  In  particular,  relaxation  functions  of  the  form 

N  -a.t 

a(t)  s-  $.e  3  ,  t  >  0  ,  (3.29) 

j-1  3 

with  >0  for  j  ■  1,2,...,N,  which  are  commonly  employed  in 

applications  of  viscoelasticity  theory,  satisfy  the  assumptions  of  Theorem 
3. 1. 

It  is  interesting  to  observe  that  if  the  relaxation  function  is  a  single 
decreasing  exponential  of  the  form  a(t)  =  e  st,  then  (2.21)  corresponds  to  a 
third  order  partial  differential  equation  without  memory.  Indeed,  in  this 
case  (2.21)  becomes 

utt(x,t)  -  ♦(ux(x,t))j£  -  «  Jq  e_a(t"T)*(ux(x,T))xdT  (3.30) 

+  g(x,t),  x  6  B,  t  >  0 
and  differentiation  of  (3.30)  with  respect  to  t  yields 

uttt(x,t)  "  ♦<Vx»t))*t  ’  «t(«x<*»t))x  (3.31) 

+  a2  Jj  e’*(t"T)*(ux(x,T))xdT 
+  gfc(x,t),  x  e  B,  t  >  0  . 

It  follows  from  (3.30)  that 

°2  Jo  e~a(t~T)*<ux<x,T))xdT  *  “♦<ux(x,t))x  (3.32) 

+  ag(x,t)  -  autt(x,t),  x  e  B,  t  >  0  . 

Substituting  (3.32)  into  (3.31)  and  using  the  definition  of  %•  obtain 

uttt  +  outt  "  ♦(ux)xt  +  aX(Vx  +  *t  +  ag  *  (3*33) 

Greenberg  [8]  studied  equation  (3.33)  with  B  ■  (0,1]  and  g  5  0  under 

homogeneous  Dirichlet  boundary  condtions.  He  derived  a  priori  estimates  which 

show  that  any  sufficiently  smooth  and  small  solution  decays  to  sero 

exponentially  as  t  ♦  •.  His  analysis  relies  on  the  Po incarl  inequality  and 
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consequently  does  not  apply  If  B  is  unbounded.  In  the  next  section,  we 
prove  Theorem  3.1  for  equation  (3.30). 

In  order  to  isolate  the  effects  of  nonlinearity  in  the  memory  term,  Hrusa 
[16]  has  studied  (2.21)  in  the  special  case  that  ♦  ia  linear  (i.e., 

♦  (C)  =  cC  for  sosw  constant  c  >  0),  but  ♦  is  allowed  to  be  nonlinear. 

His  results  apply  to  initial-boundary  value  problems  as  well  as  pure  Initial 
value  problems.  It  is  shown  in  [16]  that  the  local  behavior  of  solutions  of 
utt(x,t)  -  cuj{x(x,t)  +  J*  a'(t-T)f(ux(x,T))xdT  (3.34) 

+  g(x,t),  x  6  B,  t  >  0  , 

is  quite  similar  to  that  of  solutions  of  the  semilinear  equation 

utt  “  “xx  +  ♦<ux)  +  9  •  (3*35) 

In  particular,  a  pointwise  bound  on  ux  is  sufficient  to  continue  a  C 

solution  u  globally.  Moreover,  if  ♦'  is  bounded,  then  (3.34)  has  globally 

defined  smooth  solutions,  even  for  large  initial  data  -  independently  of  the 

sign  of  the  memory  term.  (This  requires  only  local  assumptions  on  a.)  Some 

decay  results  for  solutions  of  (3.34)  which  allow  the  data  to  be  large  are 

also  established  in  [16]. 

Several  authors  have  analysed  the  similar  first  order  problem 

ufc(x,t)  +  ♦<u(x,t))j{  +  Jq  a'  (t-T)f(u(x,T))xdt  ■  0,  (3.36) 

x  6  B,  t  >  0  , 

u(x,0 )  *  Uq ( x)  ,  x  e  B  .  (3.37) 

Equation  (3.36)  is  simpler  than  (2.21)  in  that  it  is  of  first  order,  yet  it 
retains  many  of  the  important  qualitative  features  of  (2.21).  The  chief 
motivation  for  studying  (3.36)  has  been  to  gain  insight  into  the  behavior  of 
solutions  of  (2.21). 

If  a'  vanishes  identically,  then  (3.36)  reduces  to  the  (scalar) 


conservation  law 


Ufc  +  #(u)x  ■  0 


(3.38) 


It  is  well  known  that  (3.38),  (3.37)  doea  not  generally  have  a  globally 
defined  smooth  solution,  no  matter  how  smooth  uQ  is.  Nohel  [28]  has  shown 
that  under  reasonable  conditions  on  4,  f,  and  a,  the  initial-boundary 
value  problem  (3.36),  (3.37),  with  B  *»  [0,1]  and  periodic  boundary 
conditions,  has  a  unique  global  smooth  solution  if  ug  is  sufficiently  smooth 
and  small.  (Here,  Ug  should  be  small  in  the  H2(0,1)  norm.)  Malek-Madani 
and  Nohel  [24]  have  studied  the  formation  of  singularities  in  smooth  solutions 
of  (3.36)  with  B  *  It.  Under  certain  assumptions  on  4.  ♦»  And  a  (which 
include  (3.7)  and  convexity  of  4) ,  they  give  rather  precise  conditions  on 
Ug  under  which  (3.36),  (3.37)  has  a  local  smooth  solution  for  which  first 
derivatives  became  infinite  in  finite  time. 

Relatively  little  is  known  about  weak  solutions  of  (3.22)  or  (3.36). 
Dafermos  and  Hsiao  [5]  have  established  existence  of  global  weak  solutions  (of 
class  BV)  to  systems  of  conservation  laws  with  memory  in  one  space  dimension 
for  initial  data  having  small  total  variation.  They  allow  for  very  general 
types  of  memory  terms  in  the  equations.  Their  global  results  apply  in  several 
situations  of  physical  interest  ( including  the  heat  flow  problem  mentioned  in 
Remark  3.4),  but  not  to  (3.22)  under  assumptions  which  are  appropriate  for 
viscoelastic  solids  of  the  Boltzmann  type.  (Their  procedure  does,  however, 
yield  local  (in  time)  existence  of  BV  solutions  to  (3.22)  in  this  case.) 

In  [22],  MaeCamy  studies  several  aspects  of  weak  solutions  of  equation 
(3.36)  with  ♦  3  He  also  discusses  global  existence  of  smooth  solutions 
for  ssuill  data,  and  the  formation  of  singularities  in  smooth  solutions. 
Greenberg  and  Hsiao  [9]  have  studied  the  Riemann  problem  for  a  system  which 
corresponds  to  (3.36)  with  a(t)  5  e  ,  a  >  0.  (See  also  [10].) 
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Th«  results  of  Coleman  and  Gurtln  [2]  on  mvi  propagation  (which  wsrs 
discussed  at  the  beginning  of  this  section)  hold  for  a  sore  general  class  of 
materials  with  memory.  For  these  materials,  the  displacement  u  obeys  an 
equation  of  the  form 

utt<x,t)  ■  G(u*(x,»)>  +  f ( x,t) ,  x  e  B,  t  >  0  ,  (3.39) 

where  G  is  a  smooth  (nonlinear)  functional  defined  on  a  function  space  of 
fading  memory  type,  and  for  each  x  e  B,  t  >  0, 

u^(x,s)  i-  ux(x,t-s)  V  s  >  0  ,  (3.40) 

i.e.,  u^  is  the  history  up  to  time  t  of  the  strain,  under  physically 
reasonable  assumptions  on  G,  Hrusa  [14]  has  established  global  existence 
(and  decay)  of  smooth  solutions  to  certain  history-boundary  value  problems  for 

(3.39)  with  B  ■  [0,1]  and  suitably  smooth  and  small  data.  See  also  [15] 

(and  the  references  therein)  for  a  more  complete  discussion  of  equetion 

(3.39) . 

4.  The  Cauchy  Problem  with  an  Exponential  Kernel 

In  this  section  we  sketch  the  proof  of  global  existence  of  smooth 
solutions  to  the  initial  value  problem  (3.22),  (3.23)  for  sufficiently  smooth 
and  small  data  in  the  special  case  that  the  relaxation  function  is  a 
decreasing  exponential  of  the  form  a(t)  =  e  .  We  also  discuss  the 
modifications  required  to  treat  more  general  relaxation  functions.  A  linear 
rescaling  of  time  shows  that  without  loss  of  generality  we  may  assume  a  ■  1. 
For  simplicity  we  take  g  s  0. 

In  particular,  we  consider  the  initial  value  problem 

utt(x,t)  -  ♦(uj{(x,t))x  -  Jq  e  (t  T,t(ux(x,t))xdx  , 

(4. 1 ) 

x  e  «,  t  >  o  , 
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(4.2) 


u(x,0)  -  uQ(x),  «fc(x,0)  ■  Uj(x),  x  6  K  . 

Obnm  that  tha  correaponding  equilibrium  atrasa  function  ia  given  by 

X(C)  i-  4(C)  -  4(0  »(ei  . 

Concerning  4,  4#  and  x  make  tha  aaauatptiona 

4,  4  e  c3(m),  4(0)  -  4(0)  -  o  , 

4*(0)  >  o#  4* <o)  >  o,  x'(0)  >  o  . 


(4.3) 


(4.4) 


(4.5) 


Propoaltloni  Aaaume  that  (4.4)  and  (4.5)  hold.  Than,  there  axiata  a  oonatant 


M  >  0  auch  that  for  each  uQ ,Uj  i  R  ♦  K  with 


“o e  « <•’<»> 


(4.6) 


J^,{u'(x)2  +  uj(x)2  ♦  uj"(x)2  ♦  u^x)2 

+  uj(x)2  +  u^(x)2}dx  <  II2  , 


(4.7) 


tha  initial  value  problem  (4.1),  (4.2)  haa  a  unique  eolution  u  :  R  *  [0,»)  ♦  ■ 


u  6  C2 (■  *  [0,»))  , 


(4.8) 


ut'ux'utt'utx'uxx'uttt'uttx'utxx'uxxx  9  L  <»>>  *  (4’9) 


Moreover,  aa  t  ♦  • 


♦  0  in  I.  (■)  , 

tt  tx  XX 


vvutt'utx'uxx  * 0  uniCormly  0,1  *  • 


(4.10) 

(4.11) 


Proof:  We  chooae  a  aufficiently  email  positive  number  6  and  modify  4  and 
4  (and  hence  alao  x)  amoothly  outaida  tha  interval  [-6,6]  in  auch  a  way 
that  4'  and  4'  era  conatant  outaida  [-26,26]  and 

1  <  4'(C)  <  4.  i  <  4'(C)  <  4.  X  <  X’(C>  <  X  v«e*  ,  (4.12) 

whara  4,  4,  J[,  4,  X  are  positive  conatanta  aatiafying 
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X2  -  ±  1  <  0  •  (4.13) 

(This  can  always  bs  accomplished  by  virtue  of  (4.3)  and  (4.5).)  There  is  no 
harm  in  making  this  modification  because  we  will  show  a  posteriori  that 
|ux(x,t)|  <  6  for  all  x  e  R,  t  >  0. 

Making  only  minor  changes  in  the  proof  of  Theorem  2.1  of  [7],  one  can 
establish  the  following  local  existence  result:  (4.1),  (4.2)  has  a  unique 
local  solution  u  defined  on  a  maximal  time  interval  [0,Tg) ,  T0  >  0,  with 

u  e  C2 (R  x  [0,Tg) )  ,  (4.14) 


u  ,u  ,u..  ,u  ,u  ,u... ,u.  ,u  ,u 

t  X  tt  tx  XX  ttt  ttx  txx  XXX 

Moreover,  if 


6  C(  [0,T0)i  lf( »)) 


(4.15) 


sup 

te[0,T0) 


J*  {u2  +  u2  +  u2  +  u2  +  u2  +  u 
J--  t  x  tt  tx  xx 


2 

ttt 


+  “ttx  +  “txx  +  “xxx>(x't)dx  < 


(4.16) 


then  T_ 


He  now  proceed  to  establish  a  priori  estimates  for  the  local  solution 
u  which  will  show  that  if  (4.7)  is  satisfied  with  |i  sufficiently  small  then 
(4.16)  holds.  For  this  purpose  it  is  convenient  to  introduce 


U0  :*  J^,{u^(x)2  +  u3<x>2  ♦  u'"(x)2  +  u1(x)2 


+  u'(x)2  ♦  u-(x)2}dx  , 


(4.17) 


B(t)  :» 


max  j^{u2  +  u 
«6( 0,t] 


2.2  .2  .2 
I  ♦  u^  +  u._  +  u 

X  tt  tx  XX 


+U2  .  ♦  u2  +  u2  ♦  u2  }(x,s)dx 
ttt  ttx  txx  XXX 


+  Ja  j"  t°2  +  U2  +  U2  ♦  U2  ♦  U2 

'0  J-w'  tt  fcX  M  ttt  ttx 


(4.18) 


♦uf  ♦  u2  }(x,s)dxds  ,  t  e  (0,TA)  , 

txx  XXX  O' 
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V(t)  : 


(4.19) 


-  sup  {u*  +  u*  +  u*  (x,s)  ,  t  6  tO,TQ)  . 
xes 

seto,t] 

Throughout  the  remainder  of  this  proof  we  use  T  to  denote  a  (possibly  large) 
generic  positive  constant  which  can  be  chosen  independently  of  uq,  u1#  and  Tq. 

The  reader  should  note  that  all  of  the  computations  which  follow  are  aimed  at 
establishing  an  a  priori  bound  of  the  form  (4.43). 

Differentiating  (4.1)  with  respect  to  t  and  substituting  for  the  integral 
term  from  (4.1)  (as  in  the  derivation  of  (3.33))  yields 

uttt  +  utt  “  ♦‘Vxt  +  *(Vx  ‘  (4‘20) 

The  required  estimates  will  be  obtained  by  combining  several  energy  identities  which 
we  derive  from  (4.20). 

We  first  multiply  (4.20)  by  u^t  and  integrate  over  space  and  time*  performing 
several  integrations  by  parts.  The  result  of  this  calculation  is 
\  +  ♦•(ux>ujx)(x,t>dx  +  jj  \t(x',)dxd* 

+  jZm  X(ux)utx(x,t)dx  -  X,<«jt>Utx(*'a,dX<to 

-  JZ.(j  utt  +  I  ♦'<Vutx  +  X<Ux)utx><K'0>dx  M.21) 

+  3  jJ  il  ♦-(ux)ujx(x,s)dxds 
V  t  e  [0,TQ)  . 

Next,  we  multiply  (4.20)  by  ut  and  integrate  as  above,  thereby  obtaining 
u2  +  W(ux)}(x,t)dx  ♦  J*  y  (ux)u2x(x,s)dxds 

♦  ufcutt(x,t)dx  -  u2t(x,.)dxds  (4.22) 

-  J^,{j  u*  ♦  W(«x>  ♦  ututt)(x,0)dx 

V  t  B  [  0 ,  T„ )  , 

where 
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(4.23) 


m(€)  *•  ij  x(n)dn  v  c  e  r  . 

Me  multiply  (4.22)  by  (1-e),  with  0  <  e  <  1,  and  add  the  resulting  equation 
to  (4.21).  After  rearranging  certain  terms  we  have 

j  +  2(1-e)ufcutt  +  (l-e)u*}(x,t)dx 

+  ♦' (ux)utx  +  Xfux)utx  +  (1”e>w(ux>^x*t><ix 

+  io  ^-»^eutt  +  1*'<ux>  “  ‘♦'<ux>Ju^x)*x*8>dxds 

(4.24) 

■  i"-1?  “« *  i  ♦,<V'4  *  «<a.,utK 

+  (l-e)tM(ux)  +  ufcu  +  ^  u^]}(x,0)dx 
♦  5  ^0  ♦"(ux>u^x(x,e)dxds  V  t  e  [0,TQ>  . 


He  note  that  for  each  t  e  (0,1)  the  first  integrand  on  the  left  hand  side  of 
(4.24)  is  a  positive  definite  quadratic  form  in  and  utt.  Moreover,  we  have 

+  (♦,(«x)  “  e+Mu^lu^  *  eutt  +  (4  “  e*,utx  (4.25) 

which  yields  an  obvious  lower  bound  for  the  third  integral  on  the  left  hand  side  of 
(4.24)  if  e  <  !/♦. 

The  second  integral  on  the  left  hand  side  of  (4.24)  merits  special  attention. 
Observe  that  (4.4)  and  (4.12)  imply 

IX«C)I  <  XICI,  M(C)  >  j  *C2  »(««  •  (4.26) 

Therefore,  we  have 


»  I  i An  -  *  5 


(4.27) 


A  simple  computation  reveals  that  the  right  hand  side  of  this  last  inequality  is 
positive  definite  in  u„  and  u^  for  e  sufficiently  smell,  by  virture  of 
(4.13). 
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Thus,  by  choosing  c  snail  enough  in  (4.24),  we  conclude  that 

♦  UK  +  +  ^x}(X't)dX 

♦  J  o  +  utxJ(x',)dxd* 

(4.28) 

<  T  J^{uj  ♦  u*  +  u£t  ♦  ujx}(x,0)dx 

♦  TIJq  Jl  ♦"(ux)u^x(x,s)dxde|  V  t  e  [0,Tq )  . 

We  observe  that  utt(x,0)  -  ♦  ' (u^ (x) )ujj(x)  by  (4.1)  and  (4.2),  and  since 
♦*  vanishes  outside  [-24,24]  we  have 

I  Jo  flm  ♦"<ux>Utx(x,s)dxdsl 

<  sup  l♦•(ux)utx(x,s)|  Jo  u^x(x,s)dxds  (4.29) 

xei 

sC[0,t] 

<  rv(t)E(t)  v  t  e  [o ,Tq)  . 

It  now  follows  fron  (4.28)  that 

*  u’t  ♦ 

+  Jj  ♦  ujx>(x,s)dxds  (4.30) 

<  ruQ  ♦  rv(t)Kt)  v  t  e  [o,t0>  . 

To  obtain  our  next  identity  we  nultiply  (4.20)  by  \tKX  and  integrate  as 
before,  thus  producing 
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(4.31) 


j"  (4  4*(u  )u2  -  u  u..}(x,t)dx 

2  x  xx  xx  tt 


+  J*  (x'(u  )u2  -  u  u  }(x,s)dxds 

'  0  '  -m  *  XXX  XX  tt 


1  i**  2 

“  2  l-m  U  (x,t)dx 


j"  {4  ♦  '(u  )U2  -  u  u  “  4  u2  }(x,0)dx 
2  x  xx  xx  tt  2  tx 


♦  I  Jo  JZ,  ♦"(uv>utvuL(x'a>dxd8  v  fc  e  r°,v  • 


For  each  e  >  0  we  have 


I  u  u  I  <  eu2  +  4-  u2  . 
xx  tt  xx  4e  tt 


(4.32) 


He  use  (4.12)  and  (4.32)  with  e  sufficiently  small  to  obtain  lower  bounds  for  the 
first  two  integrals  on  the  left  hand  side  of  (4.31),  and  we  majorize  the  right  hand 
side  as  before.  This  yields  the  estimate 


iZm  Uxx(X,t)dX  +  Jo  Jlo  u^x(x*s)dxds 

“  r  /:.{u2t  +  u2x)(x,t)dx  -  r  Jg  u2t(x,s)dxds  (4.33) 

<  ruQ  +  rv(t)E<t)  v  t  e  to,TQ)  . 

Combining  (4.30)  and  (4.33)  we  conclude  that 

J^,.{u2  +  u2  +  u2  +  u2  +  u2  )(x,t)dx 
1  -«•  t  x  tt  tx  XX 

+  i o  J^»^utt  +  utx  +  Uwc*<X,8>d*dS  (4.34) 

<  ro0  +  rv(t)E(t)  v  t  e  Io,t0)  . 

We  now  must  obtain  similar  bounds  for  third  order  derivatives  of  u.  In  order 
to  avoid  purely  technical  complications  and  highlight  the  main  ideas,  we  give  only 
formal  derivations  of  the  remaining  energy  identities  (4.36),  (4.37),  and  (4.39). 

The  difficulty  is  that  the  local  solution  is  not  smooth  enough  to  justify  our  formal 
procedure.  However,  all  of  these  identities  are  in  fact  valid  for  our  local 


t  i-  2 


(4.33) 


(4.34) 
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solution.  They  can  be  derived  rigorously  by  approximation.  (One  way  to  do  this  is 
to  use  difference  operators.  See  (7]  for  more  details.) 


Differentiation  of  (4.20)  with  respect  to  x  yields 
uxttt  +  uxtt  "  +(ux)xxt  +  X(Vxx  • 

We  first  multiply  (4.35)  by  uxtt  and  integrate  as  before.  The  outcome 
computation  is 

4  J  {u*  +  ♦'(«  )u2  }(x,t)dx 

2  ttx  x  txx 

+  ^0  ^  Uttx(x'8)dxdS 

+  X’(ux)uxxutxx(x't)dx 

■  i5  C  X'(ux,utxx(x'8)dxds 

=  J~{i  Uttx  +  I  ♦,(Ux,Utxx  +  X‘(ux,uxxutxx}(x'0)dx  < 

<ux 'UxxUtxxUttx  +  ^  (ux*utxuttxuxxx 

+  ♦,"<u  )U2  U  U  .  ♦  4  ♦  "(u  )u  u2 

x  xx  tx  ttx  2  x  tx  txx 

+  X"(%J«txuxxutxx}(x'8)dxd8  V  *  e  [0'V  • 

Next,  we  multiply  (4.35)  by  uxt  and  integrate  as  usual  to  obtain 


4  J* (u2  +  X'(u  )u2  } (x,t)dx 

2  -•  tx  x  xx 


+  J*  ♦,(u  )u2  ( x, s ) dxds 

'0  '-«•  x  txx 


+  Utxutxx(x't)dx  •  ->0  C.  uttx(x'8)dxds 
\x  +  i  X’(VUxx  +  utxuttx}(x'0)dx 


+  J  A  I  U"(u  )u  u  u 
'0  T  x  tx  XX  txx 


+  ?  X"(ux)utxuxx}(x'8)dxds  v  t  e  t°'To)  • 


4.35) 
of  this 


4.36) 


4.37) 
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Taking  a  suitable  linear  combination  of  (4.36)  and  (4.37)  and  estimating  the  left 
hand  side  from  below  and  the  right  hand  side  from  above  as  we  did  with  (4.21)  and 
(4.22)  shows  that 


J  (uj  +  u2  +  u2  +  u2  }(x,t)dx 
1  -m  tx  xx  ttx  txx 


+  ^  +  \xx}(x'8)dxd8 


(4.38) 


<  ruQ  +  r{v(t)  +  v(t)2}E(t)  v  t  e  [o,t0)  . 


To  obtain  our  final  identity,  we  multiply  (4.3S)  by  uxxx  and  integrate  as  usual. 
This  yields 


j"  (r-  ♦'(u  )u2  -  u  u  }(x,t) 

-*•  2  x  XXX  XXX  ttx 


dx 


♦  in  j"  (x'(u  )ux  -  U  U..  } ( X , s ) dxds 

J  0  1  — «•  v  vyy  yyv  ttx 


X  XXX  XXX  ttx 


1  .  m  o 

o  J __  u.  (x,t)dx 

2  -«•  txx 


txx 


I  (r  4 1  (u  )u^  —  u  U  *  x  u^  }(x,0)dx 

1  -•  2  W  x  xx  xxx  txx  2  txx  ’ 


1  2 

'x“xx  “xxx“txx  2  Utxx‘ 


+  in  J  Cr  ♦  "(u  )u  u2  -  24"(u  )u  u  u 

J  0  J  -m  2  *  x  tx  xxx  *  X  XX  txx  XXX 


(4.39) 


-  ♦"(u  )u.  U2  -  ♦,,,(u  )u  U2  u 

X  tx  XXX  X  tx  XX  xxx 

-  x"(u  )u2  u  } ( x , s ) dxds  V  t  e  [0 , Tn ) 

X  XX  xxx  0 

Treating  (4.39)  in  the  same  fashion  as  (4.31)  and  combining  the  result  with  (4.38) 
gives 

j"  {u2  +  u2  +  u2  +  u2  +  u2  } (x,t)dx 

J  -•  tx  XX  ttx  txx  xxx 


+  J  Jr  j"  {u2  +  u2  +  u2  }(x,s)dxds 

■'O  1  ttx  txx  xxx 


(4.40) 


<  ruQ  +  r{v(t)  +  v(t)  >E(t)  v  t  e  [o,tq)  . 
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Squaring  (4.20)  we  get 


u2  4  44 ' (u  )u2  +  44"(u  )u2  u2 

ttt  x  txx  v  x  tx  xx 


+  4uf.  +  4x'(u  )2u2  , 

tt  *  X  XX 


(4.41) 


which,  in  conjunction  with  (4.34)  and  (4.40)  yields 


JZ,  u2tt(x,t)dx  +  jJ  u2tt(x,e)dxds 

<  ruQ  +  r{v(t)  +  u(t)2}E(t)  v  t  e  [o,tq)  . 

Combining  (4.34),  (4.40),  and  (4.42)  we  finally  deduce  that 

E(t)  <  ru0  ♦  T{v(t)  +  v(t)2>E(t)  v  t  e  to,T0)  , 


(4.42) 


(4.43) 


where  T  denotes  a  fixed  positive  constant  which  can  be  chosen  independently  of 
u0,  u1 ,  and  T„. 

At  this  point,  it  should  be  noted  that  there  are  several  other  ways  to  obtain 
seme  of  the  estimates  leading  to  (4.43).  In  particular,  (4.1)  can  be  used  to 
express  uxx  in  terms  of  utt  (and  "small"  correction  terms)  through  an  inverse 
linear  Volterra  operator.  This  eliminates  the  need  for  the  identities  (4.31)  and 


(4.39). 

He  are  now  ready  to  synthesize  the  proof.  He  choose  E,  y  >  0  such  that 
_  —  1 A  —  1—  2  —  2  1  — 

r{(2E)^+  2E)  <  j,  E  <  6  ,  and  Ty  <  E.  Suppose  now  that  (4.7)  holds  with  the 
above  choice  of  y.  It  follows  from  the  Sobolev  embedding  theorem  that 

V(t)2  <  2E(t)  V  t  e  [0, TQ>  .  (4.44) 

Therefore,  we  conclude  from  (4.43)  that  for  any  t  e  [0,Tg)  with  E(t)  <  E  we 
actually  have  E(t)  <  j  E.  Consequently,  by  continuity, 

E(t)  <  jE  Vt  e  C0,Tfl )  ,  (4.45) 


provided  that 


E(0 )  <  ~  E. 
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If  necessary,  we  can  always  choose  a  smaller  u  >  0  such  that  (4.7)  implies 
E(0)  <  •j  E  and  hence  also  that  (4.45)  Is  satisfied.  This  immediately  yields 
Tp  *  •  by  virture  of  (4.16).  It  then  follows  froa  (4. It),  (4.45),  and  standard 
embedding  inequalities  that  (4.10)  and  (4.11)  hold.  Finally,  we  note  that  since 
E  <  62,  (4.19),  (4.44),  and  (4.45)  show  that  |ux(x,t)|  <  4  for  all 

x  e  R,  t  >  0.  This  completes  the  proof.  ■ 

We  close  with  a  few  remarks  concerning  the  modifications  required  to  treat  the 
Cauchy  problem  with  a  more  general  relaxation  function.  As  noted  earlier,  estimates 
for  certain  higher  order  derivatives  can  be  obtained  directly  from  (3.22)  using  the 
procedure  of  Dafermos  and  Nohel  [7].  Under  the  assumptions  of  Theorem  3.1,  equation 
(3.22)  can  be  written  in  the  form 


uttt  +  a<0)"1utt  *  ^Vxt  +  a(0)"1x(ux,x 


+  a(0 )_1f  +  tK*(4(ux)x  +  f  -  utt)l  , 


(4.46) 


where  the  *  denotes  convolution  (with  respect  to  the  time  variable)  on  [0#t] , 


i  e  6  e 

(v*w)(t)  jjj  v(t-T)w(T)dT  ,  t  >  0  ,  (4.47) 

and  K  is  the  solution  of  a  certain  integral  equation  involving  a'  and  a". 

Equation  (4.46)  is  quite  similar  to  (4.20)  and  this  suggests  a  natural 
procedure  to  get  estimates  for  the  lower  order  derivatives,  nie  chief  difficulty 
lies  in  handling  the  convolution  term.  Hiis  is  accomplished  by  rewriting  it  in 
several  convenient  equivalent  forms  involving  derivatives  on  which  we  already  have 
information.  The  details  are  carried  out  in  [17]. 
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